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I.  INTRODUCTION 


This  report  is  concerned,  with  seme  of  the  practical  aspects  of  Shanr._.kfs 
Theorem,  i.e.,  that  a hand-limited  time  series  is  completely  determined  by 
its  value  at  equally  spaced  time  intervals,  At,  if  l/2At  is  greater  than  the 

bandwidth.  Thus,  . 




f(t)  = ^ f(nAt)  h(t  - nAt)  , 


where  l/2At  > bandwidth  of  f.  5 0oe<  * \ ) 

I y 

In  practice,  f is  only  approximately  band  limited  and  the  series  contains 
only  a finite  number  of  terms.  In  this  report,  the  effects  of  bandwidth  and 
truncation  of  the  series  are  studied. 


The  programing  and  other  assistance  of  Miss  Cheryl  A.  Burr  are  hereby 
gratefully  acknowledged. 


II.  EXPLANATION  OF  SHANNON'S  THEOREM 


Shannon's  Theorem  states  that  if  x(t)  is  a time  series  with  center 

frequency  co  and  bandwidth  B,  then 
o ’ 


where 


*<t>  - X x(rn-)  h(t  - n-r)  , 


, sin  ■*- 

h(t)  = cos  g (2m  - 1) 


xt  » 
2t 


t i i- 

2B  * 


o)q  s (2m  - l)  ^ 


The  parameters  of  PROGRAM  RANDOM  that  occur  above  are  as  follows: 


T = t 


TALI  S T 
PI  « « 
M = m 


It  is  evident  from  the  form  of  the  function  f that  the  heaviest 
contributors  to  the  sum  will  be  those  for  which  nT  is  near  t.  This  observa- 
tion is  enployed  in  this  study  in  that  the  series  is  always  evaluated  only 
for  those  n such  that  t - nT  is  small.  Thus,  the  values  of  f(nT)  used  are 
those  with  nT  in  an  approximately  synmetric  interval  about  t. 
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III.  EXPLANATION  OF  PARAMETERS  READ  IN  FROM  DATA  CARD 


A.  NN  - used  to  indicate  vhich  plotting  routines  are  to  be  used.  If 
NN  = 1,  only  the  calculated  graph  will  be  plotted.  If  NN  = 0,  both 
the  original  graph  and  the  calculated  graph  will  be  plotted. 

B.  NOY  - the  "number  of  Y's",  or  number  of  data  points  used.  In  PROGRAM 
RANDOM  using  Magnetic  Tape  138,  NOY  = 3000. 

C.  LP  - the  division  factor  used  to  calculate  or  reconstruct  the  data 
curve.  When  LP  = 3,  every  third  point  of  the  data  is  used  to  recon- 
struct the  curve.  When  LP  = 5,  every  fifth  point  is  used. 

D.  ITT  - the  size  of  the  data  block  used  to  reconstruct  the  data  points. 
ITT/2  is  the  number  of  points  taken  on  each  side  of  the  reconstructed 
point  used  for  calculation.  ITT  has  been  equal  to  10,  20,  30,  hO,  100, 
and  200. 

E.  I£UM  - the  number  of  data  points  to  be  reconstructed.  LSUM  may  be  any 
integer  less  than  or  equal  to  NOY.  It  is  usually  set  at  40. 

F.  M - an  integer  for  which  the  band- limited  signal  is  between  (M  - l)  • W 
and  M • W Hz . When  LP  = 3,  M = 2.  When  LP  = 5,  M = 3. 

G.  DELTA  - the  number  of  seconds  between  each  data  point.  DELTA  = 0.00005. 

H.  NNA  - the  first  point  of  the  data  block  to  be  reconstructed.  Here 

NNA  = 1330.  When  the  data  are  unfiltered  (analog  filtered)  the  parameter 
NA  is  set  equal  to  NNA.  When  the  data  are  filtered,  NA  is  set  equal  to 
NNA  - (FPL  - l)/2  since  the  filtered  data  are  moved  over  to  the  left  by 
the  number  of  the  filter  weights  divided  by  2. 
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NNZ  - the  last  point  of  the  data  block  to  be  reconstructed.  NNZ  = 1370 
when  LP  = 5.  NNZ  = 1372  when  LP  = 3 for  purposes  of  division.  When 
the  data  are  unfiltered  (analog  filtered),  the  parameter  NZ  is  set 
equal  to  NNZ.  When  the  data  are  filtered,  NZ  is  set  equal  to  NNZ  - 
(FPL  - l)/2  since  the  filtered  data  are  moved  over  to  the  left  by  * 3 
number  of  the  filter  weights  divided  by  2. 

FPL  - the  number  of  filter  weights  used.  In  these  programs,  BPF1 
(Band-pass  Filter  1,  3 kHz-7  kHz)  and  BPF2  (Band-pass,  Filter  2, 

4 kHz-6  kHz)  both  have  201  filter  weights. 


IV.  EXPLANATION  OF  PARAMETERS  CALCULATED  IN  PROGRAM  RANDOM 


The  following  parameters  are  in  the  order  in  which  they  are  calculated 

in  PROGRAM  RANDOM. 

A.  ND  - number  of  data  points  used  from  tape. 

B.  TAU  - self-explanatory;  with  varying  data  the  parameter  TAU  equals  the 
division  factor  (LP)  multiplied  by  DELTA. 

C.  NA  - (see  Chap.  Ill,  H. ). 

D.  NZ  - (see  Chap.  Ill , I.  ) . 

E.  S Array  - Y is  an  integer  array  and  must  be  converted  to  real  and  stored 
in  the  array  S for  use  in  SUBROUTINE  FILTER. 

F.  XT  Array  - for  saving  storage  space. 

G.  N - the  number  of  weights  put  in  as  data.  These  must  have  a mirror 
image  made  in  order  to  get  the  total  array  of  weights. 

H.  NP1  - Since  there  are  N weights  read  in,  their  subscripts  will  range 
from  1 to  NP1  (or  N + l). 

I.  NFMP1  - the  total  number  in  the  array  of  weights  after  having  been 
mirror-imaged. 

u.  W Array  - the  weights. 

K.  NOX  - the  number  of  points  in  the  array  to  be  used  for  calculation  which 
equals  the  total  number  of  points  divided  by  the  division  factor  LP. 
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L.  Z Array  - the  array  consisting  of  every  third  point  (LP  = 3)  or  every 
fifth  point  (LP  = 5 ) of  the  original  data.  The  Z array  is  used  for 
calculating  the  complete  array. 

M.  T - the  time  (equals  I*DELTA) . 

N.  R - parameter  that  gives  the  subscript  of  the  element  in  the  Z array 
making  that  element  equivalent  to  the  corresponding  element  in  the 

Y array. 

O.  LSTART  - the  subscript  of  the  initial  element  used  in  the  Z array  in 
scanning  ITT  number  of  points  for  calculating  one  point.  In  this 
scanning  the  subscripts  range  from  LSTART  to  LSTOP. 

P.  LSTOP  - (see  0.). 

Q.  ZMEAJJ  - the  mean  that  is  to  be  subtracted  from  every  point  in  the  Z 
array. 

R.  TAR  - taken  directly  from  the  formula  (see  Chap.  II). 

S.  A - taken  directly  from  the  formula  (see  Chap.  II.). 

T.  B - taken  directly  from  the  formula  (see  Chap.  II.). 

U.  C - taken  directly  from  the  formula  (see  Chap.  II.). 

V.  S Array  - the  newly  calculated  array. 

W.  N1  - the  sup  error,  i.e.,  the  greatest  deviation  between  any  two 
corresponding  points  in  the  original  array  and  the  calculated  array. 


10 


X.  IWS  - the  rms  error,  i.e.,  the  average  deviation  over  all  the  points 
from  the  original  data  to  the  calculated  data. 

Y.  0 Array  - the  original  array  to  be  used  for  plotting. 

Z.  C Array  - the  calculated  array  to  be  used  for  plotting. 

AA.  X Array  - the  value  of  the  subscripts  I,  I = 1 through  ND,  to  be 
used  for  plotting  both  against  the  0 array  and  the  C array. 
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VIII.  A.  EXPLANATION  OF  GRAPHS 


The  graphs  hound  within  the  text  of  this  report  are  reproduced  data 
calculated  using  Shannon's  Theorem.  In  an  envelope  at  the  hack  of  the 
report  are  three  transparencies.  The  first  of  these  shows  the  original 
data,  which  has  heen  analog  hand-pass  filtered  (4-6  kHz)  and  digitized  at 
20  kHz.  The  graphs  of  Chapter  VIII,  B.  are  attempts  at  reproducing  these 
data.  As  is  easily  seen  from  overlaying  the  graphs  and  from  the  peak  and 
rms  errors,  these  reproductions  are  not  very  accurate.  The  second  and  third 
transparencies  and  the  graphs  of  Chapter  VIII,  C.  and  VIII,  D.  are  the  re- 
sults of  digitally  filtering  the  data  on  the  first  transparency,  using  hand- 
pass  filters  BPF1  and  BPF2  respectively,  and  then  applying  Shannon's  Theorem. 

Data  from  outside  the  displayed  interval  were  used  in  the  calculations 
when  necessary. 
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Effect  of  Shannon's  Theorem  Used  with 
Analog  Filtered  Data,  Illustrated  by  Graphs 


REPRODUCED  WAVEFORM  WITH  No.  OF  POINTS  FROM  EACH  SIDE  USED  FOR  CALCULATK'N 

ANALOG  FILTERED 


REPRODUCED  WAVEFORM  WITH  No.  OF  POINTS  FROM  EACH  SIDE  USED  FOR  CALCULATION 

ANALOG  FILTERED 
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C.  Effect  of  Shannon's  Theorem  with  Data  Filtered 
By  BPF1  (3  kHz-7  kHz),  Illustrated  by  Graphs 
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VIII. 

C.  1.  Using  Every  Third  Point  (LP  - 3) 


REPRODUCED  WAVEFORM  WITH  No.  OF  POINTS  FROM  EACH  SIDE  USED  FOR  CALCULA  f'ON 

FILTERED  BY  BPF1  (3  kHz  - 7 kHz) 


REPRODUCED  WAVEFORM  WITH  No.  OF  POINTS  FROM  EACH  SIDE  USED  FOR  CALCULATION 

FILTERED  BY  BPF1  (3  kHz  - 7 kHz) 


VIII. 

C.  2.  Using  Every  Fifth  Point  (LP  = 5) 
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VIII. 

D.  Effect  of  Shannon’s  Theorem  vith  Data  Filtered 
By  BPF2  (4  kHz-6  kHz),  Illustrated  by  Graphs 
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REPRODUCED  WAVEFORM  WITH  No.  OF  POINTS  FROM  EACH  SIDE  USED  FOR  CALCULA . ION 

FILTERED  BY  BPF2  (4  kHz  - 6 kHz} 


REPRODUCED  WAVEFORM  WITH  No.  OF  POINTS  FROM  EACH  SIDE  USED  FOR  CAlCULAIiCN 

FILTERED  BY  BPF2  (4  kHz  - 6 kHz) 
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VIII. 

D.  2.  Using  Every  Fifth  Point  (LP  = 5) 
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IX. 

COMPARISON  OF  RELATIVE  ERRORS  INVOLVING  COMBINATIONS 
OF  VARIOUS  TYPES  OF  FILTERS,  THE  NUMBER  OF  POINTS 
USED  FROM  EACH  SIDE  FOR  CALCULATION,  AND  THE  DIFFERENT 
VALUES  OF  LP,  ILLUSTRATED  BY  GRAPHS  AND  TABLES 
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RELATIVE  ERROR  WHEN  LP  EQUALS  5 
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CAB  • BJB 
6-3-67 
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TABLE  I 


RELATIVE  ERRORS 


when  LP  = 3 


FILTER 

NO. 

OF  PTSX. 

FROM  EAClK 
SIDE  USED  N. 

ANALOG 

FILTERED 

BPF1 

(3  kHz-7  kHz) 

BPF2 

(4  kHz-6  kHz) 

max  = 0 . 182 

max  = 0.088 

max  = 0.086 

3 

rms  = 0.082 

rms  = 0.054 

rms  = O.O52 

max  = 0.165 

max  = 0.059 

max  = 0.053 

5 

rms  = O.O78 

rms  = 0.034 

rms  = 0.033 

max  = O.I78 

max  = 0.034 

max  = 0.033 

10 

rms  = 0.071 

rms  = 0.020 

rms  = 0.019 

max  = 0.148 

max  = 0.025 

max  = 0.020 

15 

rms  = O.O67 

rms  = 0.011 

rms  = 0.011 

max  = 0.144 

max  = 0.021 

max  = 0.020 

20 

rms  = 0.068 

rms  = 0.013 

rms  = 0.014 

max  = 0.182 

max  = 0.013 

max  = 0.008 

100 

rms  = 0.077 

rms  = O.OO5 

rms  = 0.005 

max  = 0 . 182 

max  = 0.004 

max  = 0.004 

200 

rms  = 0.077 

rms  = 0.003 

rms  = 0.003 

max  is  the  maximum  relative  error, 
/•ms  is  the  relative  rms  error. 
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TABLE  II 


RELATIVE  ERRORS 


when  LP  = 5 


N.  FILTER 

1 

NO.  N. 

OF  PTS  N. 
FROM  EACH\ 

ANALOG 

BPF1 

BPF2 

SIDE  USED 

FILTERED 

(3  kHz-7  kHz) 

( 4 kHz-6  kHz) 

3 

max  = 0 . 182 

max  = 0.092 

max  = O.O9O 

rms  = 0.092 

rms  - 0.062 

rms  = 0.060 

5 

max  = O.186 

max  = 0.063 

max  = O.O53 

rms  = 0.077 

rms  = 0.040 

rms  = O.O39 

10 

max  = 0.148 

max  = 0.054 

max  = 0.029 

rms  = O.O69 

rms  = 0.017 

rms  = 0.015 

15 

max  = 0.153 

max  = 0.029 

max  = 0.025 

rms  = 0.070 

rms  = 0.019 

rms  = 0.016 

20 

max  = 0.148 

max  = O.O38 

max  = 0.025 

rms  = 0.070 

rms  = 0.0l4 

rms  = 0.011 

100 

max  = 0.148 

max  = 0.025 

max  = 0.012 

rms  = 0.074 

rms  = 0.009 

rms  = 0.005 

200 

max  = 0.140 

max  = 0.021 

max  = 0.008 

. 

rms  = 0.073 

rms  = 0.009 

rms  = 0.004 

max  is  the  maximum  relative  error. 


rms  is  the  relative  rms  error. 
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X.  CONCLUSION 

The  results  of  this  study  are  relatively  consistent  with  the 
anticipated  results.  If  the  data  are  assumed  to  be  in  the  4-6  kHz  band, 

LP  = 5>  the  results  are  not  as  accurate  as  if  the  data  are  assumed  to  be 
in  the  3 .3-6. 7 kHz  band,  LP  = 3-  Accordingly,  as  more  data  points  are 
used  to  compute  the  value  of  the  time  series  at  a given  point,  the  error 
tends  to  decrease. 

The  analog  filtered  data  appear  to  have  contained  a dc  bias  and 
possibly  some  other  noise  outside  the  desired  frequency  band.  This  could 
have  occurred  after  filtering,  before  digitization,  or  during  digitization. 
The  improvement  produced  by  digitally  filtering  the  data  is  marked,  and  to 
some  extent,  expected  because  of  the  low  error  level  in  digital  operations. 
The  close  correspondence  of  the  results  using  BPF1  and  BPF2  indicates  that 
most  of  the  noise  which  entered  the  digitized  data  was  outside  the  3-7  kHz 
band. 

In  summation,  there  is  a notable  correlation  between  the  original  and 
the  calculated  data  when  only  a few  data  points  are  used  in  the  calculation 
of  the  time  series.  Moreover,  there  is  a remarkable  accuracy  in  locating 
the  zeroes  of  the  functions  when  straight  line  approximations  for  both  the 
raw  data  and  the  reproduced  function  are  used. 
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